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Abstract 

We  have  computationally  designed  new  materials  for 
use  in  Naval  SOund  NAvigation  Ranging  (SONAR).  Our 
quantum-mechanical  studies  show  that  these  lead-free, 
non-toxic  oxides  will  be  high-performance  piezoelectrics 
with  promise  for  use  in  Naval  SONAR  and 
communications  applications.  To  enable  this  research, 
we  also  present  techniques  for  greatly  accelerated 
modeling  of  oxide  materials.  We  show  that  a simple 
atomistic  model  accurately  reproduces  our  quantum- 
mechanical  results  yet  is  thousands  of  times  faster.  We 
also  report  successful  porting  and  performance  timing  of 
our  computer  codes  to  the  CRA  Y XI,  resulting  in  a great 
speed-up  over  previous  architectures. 

1.  Introduction 

Perovskite  oxides  (with  formula  AB03)  have  a wide 
range  of  structural,  electrical,  and  mechanical  properties, 
making  them  vital  materials  for  many  technological 
applications,  such  as  ultrasound  machines,  cell  phones, 
and  computer  memory  devices.  Perovskite  solid  solutions 
with  a high  piezoelectric  response  are  of  particular 
interest  as  they  can  be  employed  as  sensors  in  SONAR 
devices.  When  such  a material  is  deformed  by 
underwater  sound  vibrations,  it  generates  an  electric  field 
which  can  then  be  interpreted  by  a computer  to  gain 
information,  such  as  depth  and  distance.  This  information 
is  crucial  for  the  defense  and  operation  of  Naval 
submarines  and  vessels.  Other  examples  of  perovskite 
oxide  applications  that  are  vital  in  the  challenging 
operating  conditions  of  the  military  include  non-volatile 
ferroelectric  memory  for  safe  storage  of  information,  and 
dielectric  materials  with  tunable  resonant  frequency  and 
high  signal-to-noise  ratio  for  communications  devices. 


Most  of  these  materials  are  complex  systems  with 
some  degree  of  disorder,  making  them]  challenging  to 
study  experimentally  and  theoretically.  However,  as  it  is 
their  complexity  which  gives  them  their  favorable 
properties,  highly  accurate  modeling  which  captures  the 
essential  features  of  the  disordered  structure  is  necessary 
to  explain  the  behavior  of  current  materials  and  predict 
favorable  compositions  for  new  materials.  Recently,  a 
combination  of  methodological  improvements  and  a rise 
in  computer  speed  has  made  first-principles  calculations  a 
viable  tool  for  understanding  these  complex  systems.  In 
particular,  the  density  functional  theory  (DFT) 
approach11-21  offers  a combination  of | accuracy  and 
computational  speed  that  can  reveal  the  microscopic 
structure  and  interactions  of  complex  systems13-51. 

Rational  design  of  new  materials  'with  improved 
properties  requires  simulations  that  are  accurate  and 
efficient  at  the  length  scale  of  the  processes  that  underlie 
the  desired  properties  and  that  are  performed  at  typical 
device  operation  temperatures.  This  creates  a “length 
scale  gap”  and  a “temperature  gap”  between  the  desired 
processes  and  the  DFT  modeling  of  ferroelectric  oxides 
that  is  typically  done  for  relatively  small  unit  cells  and  at 
low  temperatures.  We  bridge  these  gaps  with  two 
strategies.  First,  we  develop  correlations  between  low 
temperature  local  properties  computed  with  DFT 
calculations  and  macroscopic  observables  such  as  the 
ferroelectric  to  paraelectric  transition  temperature.  We 
also  develop  an  atomistic  model  based  on  first-principles 
calculations  that  allows  us  to  do  accurate  simulations  at 
finite  temperatures  and  at  the  nanoscale.  Using  the 
discovered  correlations  between  microscopic  and 
macroscopic  oxide  properties,  we  computationally  design 
new  environmentally  friendly  high  performance 
piezoelectric  and  dielectric  materials.  We  also  report  our 
success  in  porting  and  performance  tuning  our  codes  for 
the  Cray  XI  platform. 
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2.  Methodology 

For  all  the  DFT  calculations  presented  here  we  use  a 
plane-wave  basis  set.  Such  a basis  set  is  complete  and 
offers  the  advantage  of  carrying  out  operations  in  both 
real  space  and  reciprocal  space  through  the  use  of  fast 
Fourier  transforms.  Calculations  are  done  using  the 
standard  LDA[6'  exchange-correlation  functional  using 
our  in-house  plane  wave  code.  To  reduce  the 
computational  cost  of  the  calculations  we  use  designed 
non-local[71,  optimized18'  pseudopotentials  to  represent  the 
interactions  of  the  nucleus  and  the  core  electrons  with  the 
valence  electrons.  Minimization  of  the  energy  with 
respect  to  the  electronic  degrees  of  freedom  is  done  using 
the  blocked-Davidson'9'  iterative  diagonalization 
procedure'101  with  Pulay  density  mixing'11'.  Ionic 
minimization  is  performed  using  a quasi-Newton 
algorithm'12'.  DFT  calculations  are  performed  on 
supercells  of  up  60  atoms,  with  a variety  of  atomic 
confgurations  examined  to  ensure  accurate  modeling  of 
the  disordered  perovskite  structures.  Molecular  dynamics 
calculations  were  performed  using  the  MOLDY 
package'13'. 


Figure  1.  Pair  distribution  function  of  DFT  simulated 
structure  of  PbSc^Wi/jOj  (solid  line!  in  comparison 
with  the  experimental  data  (dots)'™1.  The  close 
agreement  demonstrates  that  our  calculations  capture 
the  behavior  of  this  complex  material.  Similar 
agreement  is  obtained  for  other  compositions  as  well. 

3.  Results 

3.1.  Structure  Property  Correlations  in  Pb  based 
solid  Solutions. 

First-principles  calculations  have  proved  to  be  widely 
useful  in  obtaining  a microscopic  understanding  of  the 
mechanisms  underlying  cation  ordering'14',  phase 
transitions'3,151  and  piezoelectricity'4,16,17'  in  complex 
ferroelectric  solid  solutions.  However,  reliable  prediction 
of  ferroelectric  to  paraelectric  transition  temperature  7"c 
from  first  principles  remains  to  be  achieved  and  the 


microscopic  mechanisms  underlying  high  Tc  (or  the 
related  temperature  maximum  of  the  dielectric  constant) 
are  currently  unclear.  Here  we  present  our  study  of  the 
correlations  between  local  structures  at  the  atomic  level  as 
determined  by  DFT  calculations  and  experimentally 
observed  transition  temperatures'18"20'.  We  perform  our 
simulations  on  PbZri^Tix03  (PZT), 

(l-x)PbMg1/3Nb2/303-xPbTi03  (PMN-PT)  and 
( 1 ~x)PbZn  i/3Nb2/303-xPbTi03  (PZN-PT)  solid  solutions 
that  form  the  core  of  current  and  next-generation  Navy 
SONAR  devices.  We  also  study  the  related  systems 
( 1 -x)PbSc2/3  W1/303-xPbTi03  (PSW-PT), 

(l-x)PbSc2/3W1/303-xPbZr03  (PSW-PZ),  and 
PbSc,/2Nb1/203  (PSN).  Our  modeling  of  local  structure  is 
accurate,  as  indicated  by  excellent  agreement  between 
pair-distribution  functions  obtained  from  DFT 
calculations  and  by  neutron  scattering'18'  (Figure  1). 

We  find  that  there  is  a strong  relationship  between 
the  low-temperature  polarization  P (as  found  by  our  DFT 
calculations)  and  the  experimentally  obtained  temperature 
Tc  at  which  the  dielectric  constant  is  maximum: 

T^yP2,  (1) 

where  y is  a constant.  Figure  2 shows  the  correlation 
between  the  transition  temperatures  predicted  by  Eq.  1 
and  those  from  experimental  data.  The  slope  of  the  least 
square  fitted  curve  guiding  the  eye  is  1.1,  indicating  that 
Eq.  1 accurately  and  quantitatively  captures  the  trends  in 
transition  temperatures  among  these  fairly  different 
systems. 

The  correlation  presented  in  Figure  2 can  be 
explained  based  on  the  simple  Landau  theory  for 
ferroelectric  crystal.  The  free  energy  is 

G = G0  + a(T-Tc)P2  + 8P*  (2) 

where  G0  is  the  free  energy  of  the  high-symmetry 
paraelectric  phase,  a and  3 are  constants  and  P is 
polarization.  The  equilibrium  polarization  at  a given 
temperature  T is 

Pz=a{Tc-T)l28.  (3) 

Substituting  T = 0 used  in  our  DFT  calculations  and 
solving  for  Tc,  we  obtain 

Tc  = IBP2/ a,  (4) 

relating  the  Tc  to  the  square  of  computed  polarization. 
The  high  quality  of  the  fit  in  Figure  2 means  that  for  the 
Pb-based  perovskites  examined  here  the  ratio  81a  is  nearly 
constant  and  can  be  approximated  by  8la= 635  Km4/C2. 
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Figure  2.  Correlation  between  transition  temperatures 
predicted  by  Eq.  1 (y  =1270  Km4/C2)  and  obtained 

experimentally  for  a variety  of  Pb-based  solid 
solutions.  The  close  match  means  that  we  can  now 
predict  transition  temperature  from  low-temperature 
calculations. 

The  linkage  between  the  microscopic  properties 
computed  with  accurate  first-principles  calculations  and 
Tc  makes  possible  the  rational  design  of  high  Tc 
piezoelectrics.  Fundamentally,  we  have  shown  that  the 
behavior  of  individual  ions  in  heterovalent  perovskite 
ferroelectric  solid  solutions  can  be  directly  linked  to  local 
structure  and  macroscopic  collective  properties 

3.2.  Bond  Valence  Model  - Efficiently  Modeling 
Complex  Oxides. 

While  some  properties  of  ferroelectric  materials  are 
due  to  the  local,  Angstrom-scale  structure,  others  are 
dependent  on  behavior  at  longer  length  and  time  scales. 
Domain  wall  motion,  crucial  for  non-volatile  ferroelectric 
memory  applications,  is  an  excellent  example  of  a 
technologically  important  but  poorly  understood 
nanoscale  process.  Highly  accurate  simulations  are 
needed  to  obtain  a quantitative  microscopic  understanding 
of  such  processes.  The  prohibitive  cost  of  a first- 
principles  calculation  on  such  a large  system  necessitates 
development  of  computationally  inexpensive,  yet  accurate 
interatomic  potentials  for  molecular  dynamics 
simulations. 

In  this  paper,  we  present  an  atomistic  model  for  the 
molecular  dynamics  study  parameterized  from  first- 
principles  calculations  and  show  the  ferroelectric  phase 
transition  of  bulk  PbTi03  crystal  and  the  domain  wall 
energy  of  180°  PbTi03  domains.  Our  potential  model 
follows  the  bond-valence  theory'211  developed  from 
experimental  crystallographic  data  for  many  oxide 
systems.  The  valence  of  an  atom  is  assumed  to  be 
distributed  among  the  bonds  it  forms,  and  the  resulting 
valence  (the  “bond  valence”)  correlates  with  the  bond 
length  through  an  inverse  power  relation.  In  other  words, 
the  sum  of  bond  valences  at  each  atom  is  equal  to  the 
atomic  valence.  The  validity  of  the  bond-valence  (BV) 


theory  for  ferroelectric  oxides  was  proved  in  our  previous 
work'3221.  . | 

The  reference  structures  were  obtained  from  the 
ab  initio  density  functional  theory  molecular-dynamics 
calculations  on  a small  supercell.  About  2000  reference 
structures  were  obtained,  and  25  of  those  were  selected 
for  the  fitting  database.  To  obtain  the  bond-valence 
model  potential  parameters,  we  used!  the  simulated 
annealing  global  optimization  method  to  minimize  the 
difference  between  the  BV  and  DFT  energies  and  atomic 
forces  for  the  structures  in  the  database,  j The  BV  model 
with  optimized  parameters  was  then  tested  against  all 
2000  reference  structures,  showing  good  j agreement  with 
DFT  results  as  shown  in  Figure  3.  | 


computed  by  DFT  (dotted)  and  by  bondj-valence  model 
calculation  with  optimized  parameters  (solid  line) 

i 

We  then  performed  molecular  dynamics  simulations 
with  the  BV  model  at  a variety  of  temperatures  in  the 
NVT  ensemble  using  the  experimental  lattice 
constants'231.  In  agreement  with  this  experiment  we  find 
that  as  temperature  increases,  polarization  P decreases. 
Using  a relatively  small  3*3*3  MD  cell,  we  find  a 
ferroelectric  to  paralectric  phase  transition  at  around 
600K.  The  phase  transition  is  of  mixed  order-disorder 
and  displacive  type,  in  agreement  with  experimental 
data'241.  As  the  size  of  the  simulation  cell  is  increased  to 
5*5*5  and  7*7*7,  the  ferroelectric-to-paraelectric  phase 
transition  temperature  increases  to  650K  and  680K 
respectively.  The  transition  temperature  for  the  largest 
simulation  cell  agrees  well  with  the  experimental  value  of 
765K.  Comparison  of  temperature  jdependence  of 
polarization  with  experimental  data  also  shows  good 
agreement'251  (Figure  4). 
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Figure  4.  Temperature  dependence  of  PbTiOj 
polarization.  Our  model  predicts  a ferroelectric  to 
paraelectric  phase  transition  at  Tc  =680K,  in  good 
agreement  with  experimental  c =765K.  Experimental 
lattice  constants  for  a and  c axes  were  used  at  each 
temperature.  Open  circles  are  the  polarization  for  the 
3x3x3  MD  cell  and  solid  triangles  are  the  polarization 
for  the  7x7x7  MD  cell. 

To  test  the  accuracy  of  our  potential  for  the  study  of 
domain  wall  motion  in  non-volatile  memory  devices,  we 
computed  the  energy  of  the  180°  domain  wall.  We  find 
that  domain  wall  energy  in  excellent  agreement  with  DFT 
data  (Table  1).  Since  our  BV  parameterization  database 
did  not  contain  any  domain  wall  structures,  the  excellent 
agreement  between  BV  and  DFT  methods  is  another 
strong  indicator  that  our  BV  model  correctly  captures  the 
essential  physics  of  PbTi03. 

Table  1.  Comparison  of  experimental,  DFT  and  BV 
data  for  PbTiOj.  Polarization  P in  C/m2,  domain  wall 
energy  in  mJ/mzand  transition  temperature  Tc  in 
Kelvin. 
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Our  development  of  a highly  efficient  yet  very 
accurate  atomistic  model  for  PbTi03  creates  an 
opportunity  for  study  of  a variety  of  technologically 
important  nanoscale  processes.  We  are  now  working  on 
techniques  to  enable  simulations  at  the  mesoscale. 

3.3.  DFT  Study  of  Non-toxic  Piezoelectrics. 

Our  previous  work  showed  that  Ag  ions  display  off- 
centering  behavior  crucial  to  high  piezoelectricity  when 
alloying  with  PbTi03  eliminates  B-06  octahedral 
rotations1261.  Building  on  our  previous  computational 
materials  design  studies,  we  further  investigated  the 
properties  of  xAgNb03-(  1 - x)PbTi03  (ANPT) 

ferroelectric  solid  solutions  as  well  as  the  properties  of  the 


lead  free,  environmentally  friendly  xAgNb03- 
(l-x)BaTi03  (ANBT)  solid  solutions. 

We  look  for  compositions  which  exhibit 
improvement  over  the  current  PZT  or  BaTi03  (BT) 
piezoelectrics  in  three  desired  properties:  piezoelectric 
performance,  Tc  and  lower  toxicity.  In  order  to  determine 
if  the  proposed  materials  are  likely  to  be  good 
piezoelectrics,  we  examine  P values  and  the  energy 
difference  between  rhombohedral  (R)  and  tetragonal  (T) 
phases  (A£R_T).  Materials  with  large  P will  couple 
strongly  to  an  applied  electric  field;  a small  A£R-T 
indicates  the  presence  of  a monoclinic  (M)  phase  at  the 
morphotropic  phase  boundary  (MPB)  crucial  for  easy 
polarization  rotation  and  high  piezoelectric 
performance14,271.  High  Tc  is  favored  by  a large  energy 
difference  between  the  distorted  polarized  structure  and 
the  non-polar  high  symmetry  structure  (ferroelectric 
instability  A£fe)  and  high  polarization,  and  high  Tc  is 
disfavored  by  large  octahedral  rotations128'291.  Elimination 
of  Pb  atoms  from  the  A-site  makes  a material 
environmentally-friendly.  The  results  for  A£R_T  are 
presented  in  Figure  5,  and  P magnitudes,  A£fh  values  and 
average  06  rotation  angles  at  MPB  of  each  solution  are 
presented  in  Table  2. 


Figure  5.  Energy  difference  between  rhombohedral 
and  tetragonal  phases  as  a function  of  AgNbOj 
content  for  ANPT  (solid)  and  ANBT  (dotted)  solid 
solutions.  Positive  AEr-t  values  indicate  that 
tetragonal  phase  is  preferred.  Morphotropic  phase 
boundaries  are  present  when  AEr.t  crosses  0 at 
37.5%  AN  content  for  ANPT  system  and  at  12.5%  and 
50%  for  ANBT  system. 
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Table  2.  Results  of  DFT  calculations  for  compositions 
of  ANPT  and  ANBT  solid  solutions  closest  to  MPB. 

Results  for  BT  and  50/50  PZT  are  given  for 
comparison.  The  new  ANBT  materials  are  lead-free 
with  piezoelectricity  on  par  with  PZT.  MPB  content  is 
in  percent,  P and  ea  magnitudes  are  in  C/m2,  AEfe 
values  are  in  meV  per  5-atom  cell,  average  0«  rotation 
angles  are  in  degrees.  The  values  for  50/50  PZT  are 
taken  from  References  17,  20,  and  26. 
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We  find  that  while  the  ANPT  system  should  exhibit 
good  piezoelectric  performance,  the  Tc  most  likely  will 
not  be  higher  than  that  of  current  PZT  materials.  Our 
calculations  predict  a standard  phase  diagram1301  with  a 
T — > M — > R sequence  of  phase  transitions  with 
decreasing  PT  content  and  a M phase  around  37.5%  AN 
composition.  The  large  0.6  C/m2  P at  MPB  is  promising 
for  good  piezoelectric  performance.  However,  Tc  at  MPB 
will  most  likely  be  significantly  the  PT  Tc  due  to  a 
decrease  in  A£fe  and  large  06  rotations  (7.7°  versus  2.0° 
in  50/50  PZT1221).  Finally,  the  presence  of  the  Pb  atoms 
on  the  A-site  makes  this  proposed  material 
environmentally  unfriendly  and  difficult  to  synthesize  in  a 
pure  bulk  form  due  to  formation  of  the  competing 
pyrochlore  phase  commonly  found  in  perovskites  that 
contain  both  Pb  and  Nb|  U|. 

Our  results  for  the  ANBT  system  indicate  that  it  is  a 
promising  candidate  for  high  performance  lead-free 
piezoelectrics.  This  system  has  a rich  phase  diagram  with 
an  unusual  R— » M — » T — » M R phase  transition 
sequence,  which  is  due  to  volume  expansion  effects.  In  a 
solid  solution,  alloying  a perovskite  of  a small  volume 
into  a perovskite  of  large  volume  expands  the  A site  for 
the  smaller  A-cation.  This  can  be  thought  of  as  applying 
negative  pressure  to  the  smaller  perovskite.  Tetragonality 
of  the  solid  solution  is  favored  both  by  a larger  negative 
pressure132,331,  (i.e.,  a larger  average  system  volume),  and 
by  a greater  fraction  of  the  smaller  A-cation  sites  under 
negative  pressure.  Since  these  two  effects  have  opposite 
dependence  on  AN  concentration,  A£R_T  can  be  expected 
to  exhibit  non-monotonic  dependence  on  AN  content,  as 
is  observed  in  our  DFT  calculations.  The  two  MPBs  are 
located  at  12.5%  and  50%  AN  content,  with  P values 
enhanced  relative  to  BT  and  comparable  to  the  P of  some 
of  the  current  single  crystal  lead  based  materials1341. 

For  both  MPB  compositions,  A£fe  values  are 
significantly  increased  relative  to  BT,  with  small  06 
rotations  for  the  12.5%  AN  composition;  for  50%  AN 


composition  the  5°  rotations  are  larger  than  2°  rotations 
found  in  50/50  PZT1221.  For  these  two  promising 
compositions,  we  also  computed  the  e33  piezoelectric 
coefficients.  The  computed  e33  values  were  9.5  C/m2  and 
14.5  C/m2  for  the  12.5%  AN  and  50%  AN  compositions 
respectively.  These  values  are  comparable  to  the  12.6 
C/m2  computed  and  11.6  C/m2  measured  for  50/50  PZT 
solid  solution1171.  In  summary,  we  have]  shown  that  Ag 
cations  can  be  ferroelectrically  active  on  the  perovskite  A- 
site.  Alloying  AgNb03  with  PbTi03  leads  to  a standard 
phase  diagram  with  an  MPB  around  37.5%  AgNb03 
composition.  An  unusual  phase  diagram'  is  obtained  for 
the  AgNb03-BaTi03  solution,  with  excellent  piezoelectric 
properties  at  the  two  MPB  locations.  The  ANBT  solid 
solutions  should  also  display  a higher  Tc  than  BaTi03. 
The  low  toxicity  and  excellent  piezoelectric  properties 
make  ANBT  solid  solutions  promising  materials  for 
experimental  investigation. 

i 

3.4.  DFT  Based  Design  of  New  Dielectric 
Materials.  | 

I 

Besides  applications  utilizing  their  ferroelectric 
properties,  perovskite  materials  are  also!  widely  use  as 
dielectrics  in  telecommunications  applications.  There  is  a 
current  need  for  materials  that  have  a 'high  tunability 
(ability  to  change  their  dielectric  constant  and  resonant 
frequency  with  an  applied  electric  field),  but  also  low 
dielectric  loss  (increased  signal-to  noise  ratio).  The 
current  state  of  the  art  tunable,  low-loss  dielectrics  are 
based  on  incipient  ferroelectrics.  These  materials  have  a 
very  shallow  ferroelectric  instability,  typically  develop 
polarization  only  at  very  low  temperature  and  are 
paraelectric  at  room  temperature.  The  flat  potential 
energy  surface  gives  these  materials  the  ability  to  develop 
a large  polarization  with  an  applied  electric  field,  leading 
to  a large  dielectric  constant.  Current  dielectrics  such  as 
KTai-xNbx03,  AgTai_xNb^03  and  Bai-jSr,Ti03  are  all 
B-site  based,  with  the  B-cations  providing  the  off- 
centering  and  polarization  under  applied  electric  fields. 
Here,  we  examine  the  possibility  of  alloying  end-member 
materials  that  are  not  incipient  ferroelectrics  to  give  rise 
to  an  A-site  based  incipient  ferroelectric  solid  solution, 
using  the  xAgNb03-(l-x)BaZr03  (ANBZ)  system  as  a 
test  case.  1 

BaZr03  is  a low-loss  dielectric  with  a volume  that  is 
22%  larger  than  that  of  AgNb03.  Recent  theoretical 
calculations  have  shown  that  large  perovskite  volume 
expansion  under  “negative  pressure”!  leads  to  an 
anomalous  increase  in  tetragonality  and  very  large  cation 
displacements132331.  Alloying  AN  with  BZ  is  a way  of 
realizing  such  negative  pressure  experimentally.  Our 
calculations  on  the  AN-BT  system  also  suggest  that 
volume  expansion  effects  are  favorable  for  ferroelectricity 
in  the  Ag-based  systems. 
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Data  in  Table  3 show  that  alloying  just  12.5%AN 
into  BZ  leads  to  a development  of  small  polarization, 
typical  of  incipient  ferroelectrics  at  very  low 
temperatures1351.  The  polarization  is  due  to  a combination 
of  large  Ag  off-centering  and  small  displacements  by  the 
other  cations.  The  large  Ag  displacements  are  due  to  a 
volume  expansion  of  the  AgOi2  cage  which  creates  a deep 
instability  in  the  Ag  potential  energy  surface.  The  small, 
but  non-zero  distortions  of  the  Ba,  Zr  and  Nb  cations 
indicate  that  the  volume-expansion  induced  Ag 
displacements  soften  the  local  potential  felt  by  the  Ba,  Zr 
and  Nb  cations,  creating  a very  shallow  double  well 
similar  to  the  shallow  double  well  found  in  SrTi03 
incipient  ferroelectric. 

Table  3.  Results  of  DFT  calculations  for  ANBZ  solid 
solution.  Displacements  in  A,  P magnitudes  are  in 
C/m2,  AEfe  values  are  in  meV  per  5-atom  cell,  average 
Os  rotation  angles  are  in  degrees. 


Ag 

disp 

Ba 

disp 

Zr/Nb 

disp 

Os  rot 
angle 

P 

AEfe 

0.0  AN 

0.00 

0.00 

0.0 

0.0 

0 

12.5  AN 

0.89 

0.06 

0.05 

3.3 

0.11 

59 

25  AN 

0.84 

0.06 

0.07 

3.4 

0.19 

64 

37.5  AN 

0.64 

0.06 

0.08 

5.9 

0.24 

34 

50  AN 

0.44 

0.07 

0.06 

8.1 

0.20 

24 

Examination 

of  the 

data 

for  the 

other 

AN 

compositions  shows  that  the  AN-BZ  solid  solution 
displays  an  interesting  C -»  T ->  M O sequence  of  0 K 
compositional  phase  transitions  and  a nonmonotonic 
dependence  of  A on  AN  content,  due  to  volume  effects 
as  explained  above  for  the  AN-BT  solid  solution.  The 
octahedral  tilt  angle  shows  an  increase  from  0°  at  pure  BZ 
to  8.1°  at  50/50  ANBZ  composition.  Large  tilt  angles  in 
perovskites  have  been  correlated  with  large  dielectric 
loss1361,  indicating  that  dielectric  loss  will  increase  with 
AN  content  as  well.  Taken  together,  the  correlations 
between  incipient  ferroelectricity  and  high  tunabilitv,  and 
between  low  dielectric  loss  and  low  octahedral  tilt  angle 
suggest  that  12.5%  and  25%  AN-BZ  compositions  are 
good  candidate  low-loss,  frequency-agile  materials. 

3.5.  Cray  XI  Port. 

During  the  course  of  the  past  year,  we  have  ported 
our  plane- wave  code  to  the  XI  platform.  Our  program 
contains  approximately  9,500  lines  of  code  and  makes 
calls  to  several  mathematical  libraries,  making  the  porting 
and  optimizing  performance  a non-trivial  task.  The 
results  for  several  benchmarks  were  tested  and  compared 
against  typical  production  platforms  such  as  Pentium  4 
Xeon  Beowulf  Cluster  and  the  SGI  03K  computers. 
Agreement  to  six  decimal  places  between  the  XI  and 


other  platforms  was  achieved.  As  expected,  proper 
vectorization  and  optimal  memory  usage  were  crucial  for 
taking  advantage  of  the  Cray  XI  architecture.  After  the 
code  was  fully  vectorized,  a speed-up  of  9.2  times  over 
Pentium  4 Xeon  and  6.6  times  over  the  SGI  03K  was 
achieved.  We  are  currently  using  the  XI  in  production 
mode  and  working  on  optimizing  scratch  disk  and 
memory  usage  to  get  further  speed-up. 
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4.  Significance  to  DoD 

The  perovskite  oxides  are  used  extensively  in  modem 
Naval  SONAR  devices,  non-volatile  memories  and 
telecommunications  applications.  The  US  Navy  would 
reap  a considerable  military  advantage  from  developing 
SONAR-detecting  materials  and  microwave  materials 
with  higher  performance,  lower  cost  and  less  harmful 
environmental  side  effects.  Understanding  the  behavior 
of  current  perovskite  oxides  is  critical  for  the  goal  of 
developing  new  materials.  Once  the  relationship  between 
the  atomic  composition,  structure  and  materials  properties 
are  understood,  new  materials  that  improve  upon  existing 
technology  can  be  designed.  Our  DFT  calculations  have 
revealed  the  microscopic  origin  of  transition  temperature 
which  sets  the  operating  range  of  piezoelectric  devices, 
and  our  development  of  the  bond-valence  modeling 
approach  will  be  of  great  benefit  to  the  study  of  processes 
on  nano-and  mesoscales  in  non-volatile  ferroelectric 
memory  and  relaxor  piezoelectrics  used  in  the  devices  of 
the  next  decade.  Our  ANBT  and  ANBZ  work  identifies 
promising  candidates  for  environmentally-friendly  high 
performance  piezoelectrics  and  dielectrics. 
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5.  Systems  Used 

ERDC  SGI3900  and  Cray  XI,  ARSC  Cray  SX-6  and 
XI,  and  AHPCRC  Cray  XI. 

6.  CTA 

Computational  Materials  Science  (CMS) 
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